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Abstract 

We develop a grid method for multi-channel scattering of atoms in a waveguide with harmonic 
confinement. This approach is employed to extensively analyze the transverse excitations and 
deexcitations as well as resonant scattering processes. Collisions of identical bosonic and fermionic 
as well as distinguishable atoms in harmonic traps with a single frequency uj permitting the center- 
of-mass (cm.) separation are explored in depth. In the zero-energy limit and single mode regime 
we reproduce the well-known confinement-induced resonances (CIRs) for bosonic, fermionic and 
heteronuclear collisions. In case of the multi-mode regime up to four open transverse channels 
are considered. Previously obtained analytical results are extended significantly here. Series of 
Feshbach resonances in the transmission behaviour are identified and analyzed. The behaviour 
of the transmission with varying energy and scattering lengths is discussed in detail. The dual 
CIR leading to a complete quantum suppression of atomic scattering is revealed in multi-channel 
scattering processes. Possible applications include, e.g., cold and ultracold atom-atom collisions in 
atomic waveguides and electron-impurity scattering in quantum wires. 
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I. INTRODUCTION 



During the last years, the field of ultracold few-body confined systems has progressed 
remarkably. By employing optical dipole traps P and atom chips [21 El S] it is possible 
to fabricate mesoscopic structures in which the atoms are freezed to occupy a single or a 
few lowest quantum states of a confining potential such that in one or more dimensions the 
characteristic length possesses the order of the atomic deBroglie wavelength. These config- 
urations can be well described by effective one-dimensional systems. Well-known examples 
are quantum wires and atom waveguides or quasi two-dimensional systems such as 2D elec- 
tronic gases. The quantum dynamics of such systems is strongly influenced by the geometry 
of the conflnement. 

To control the coherent propagation of particle beams through ID or 2D traps it is crucial 
to very well understand the impact of the conflnement on the coUisional properties (see refs. 

[siEiiTisiaiToiiiiiiiaiiaiiiiiiansiiniiTHiiB 

refs. therein). Free-space scattering theory is no longer valid in such systems and a new theory 
is needed. This stimulated the development of quantum scattering theory in low dimensions. 
For a sufficiently dilute gas under a strong transverse confinement, one may expect both the 
chemical potential and the thermal energy ksT being less than the transverse level spacing. 
With this assumption, the dynamics of ultracold atoms in low dimensional structures e.g. 
tight wave guides has been studied using the simplification that the atoms occupy only 
the ground state of the transverse confining potential. Nevertheless, the virtual transverse 
excitations in the course of the collision process can play a crucial role in the scattering 
leading to the so-called confinement-induced resonances (CIRs), predicted by Olshanii [131 
[H]. S-wave zero-energy scattering of bosons is mapped to an effective longitudinal zero- 
range pseudopotential gi^Slz) approximating the ID atom-atom interaction in a transverse 
harmonic confinement. It was shown [131 E] that CIRs appear if the binding energy of 
the pseudopotential, approximating the two-atom molecular state in the presence of the 
confinement, coincides with the energy spacing between the levels of the transverse harmonic 
potential. In the vicinity of the CIR the coupling constant gio can be tuned from —oo 
to +00 by varying the strength of the confining potential over a small range. This can 
result in a total atom-atom refiection, thereby creating a gas of impenetrable bosons. CIRs 
have also been studied for the three-body [T^l [IS] and the four-body p7| scattering under 
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confining potential, as well as for a pure p-wave scattering of fermions [18]. Experimental 
evidence for the CIRs for bosons [121 [20] and fermions [21] has recently been reported. A 
general analytical treatment of low-energy scattering under action of a general cylindrical 
confinement involving all partial waves and their coupling, was first provided in ref . |22j 
for a spherically symmetric short-range potential. The effect of the cm. motion on the 
s-wave collision of two distinguishable atoms (i.e. a two-species mixture) in a harmonic 
confinement as well as for two identical atoms in a non-parabolic confining potential has 
been investigated in ref. [23] in the zero-energy limit neglecting the s and p wave mixing. A 
detailed study including the effect of the cm. nonseparability and taking into account the 
s and p wave mixing for harmonic confinement was performed in the single-mode regime in 
ref. [26]. Recently a so-called dual-CIR was discovered [211 EHl ES], which is characterized 
by a complete transmission (suppression of quantum scattering) in the waveguide due to 
destructive interference of s and p waves although the corresponding collisions in free space 
involve strong interactions. 

The problem of atomic pair collisions under the action of a harmonic trap in the multi- 
mode regime when the energy of the atoms exceeds the level spacing of the transverse 
trapping potential is much more intricate than the single-mode regime due to several open 
transverse channels. It demands the development of a multi-channel scattering theory ac- 
counting for the possible transitions between the levels of the confining potential. Using 
as a starting point the formalism for scattering in restricted geometries suggested in refs. 
[23 I2H] , the multi-channel scattering problem for bosons in a harmonic confining potential 
has been analyzed analytically by Olshanii el al. |i29j in the s-wave pseudopotential ap- 
proximation and the zero-energy limit. Without detailization of the interatomic interaction 
but using only two input parameters - the s-wave two-body scattering length in free space 
and the trap frequency - they have derived an approximate formula for the scattering ampli- 
tude describing two boson collisions confined by transverse harmonic trap in the multi-mode 
regime. 

In the present work we develop a general grid method for multi-channel scattering of 
identical as well as distinguishable atoms confined by a transverse harmonic trap. The 
method applies to arbitrary atomic interactions, permitting a rich spectral structure where 
several different partial waves are participating in the scattering process or even the case of 
an anisotropy of the interaction. The only limitation is that we consider harmonic traps with 
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a single frequency for every atom causing a separation of the cm. and relative motion. With 
our approach we analyze transverse excitations/deexcitations in the course of the coUisional 
process (distinguishable or identical atoms) including all important partial waves and their 
couplings due to the broken spherical symmetry. Special attention is paid to the analysis of 
the CIRs in the multi-mode regimes for non-zero collision energies, i.e., to suggest a non- 
trivial extension of the CIRs theory developed so far only for the single-mode regime and 
zero-energy limit. 

In detail we proceed as follows. Section II contains the derivation of the Hamiltonian, 
definitions of the interatomic interaction and the scattering asymptotics in the confined 
geometry, and a discussion of the scattering parameters and the transition probabilities 
characterizing the two-body collisions in the trap. Our computational method is outlined 
in Sec. Ill, important technical details are given in the Appendix. In section IV our results 
are presented and analyzed. A summary and conclusions are given in Sec. V. 



II. HAMILTONIAN AND TWO-BODY SCATTERING PROBLEM IN A WAVEG- 
UIDE 

Let us consider collisions of two atoms under the action of the transverse harmonic 
confinement. We address both cases, distinguishable and indistinguishable atoms under the 
action of the same confining potential, i.e. the trap is characterized by a single frequency cu 
for every atom. The corresponding Hamiltonian is given by 

fi^ h"^ 1 1 

where nii is the mass of the i — th atom, V{ri — is the two-body potential describing the 
interaction between two colliding atoms in free space and Tj = (pj, Zi) — (r^, 9i, (pi) are the 
coordinates of the i — th atom. The Hamiltonian is separable with respect to the relative 
r = ri — r2 and cm. R = (miri -|- m2r2) / (mi -|- 777.2) variables 

H — HcM + Hjiei , (2) 

where 

Hcm^-^VI+^Mu'pI, (3) 
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and 

Hnei=-^Vl + ll^u;'p' + V{r) . (4) 

Here M = mi + m2 and fi = 17111712/ {mi + 1712) are the total and reduced masses respectively. 
The problem, thus, reduces to scattering of a single effective particle with the reduced mass 
/i and collision energy e > fkjj, off a scatterer V{r) at the origin, under transverse harmonic 
confinement with frequency uj 

|v,2 + ^/xa;V + \/(r) ^{v) = ei,{v). (5) 

Here the energy of the relative two-body motion e = e_L + ey is a sum of the transverse e±_ 
and longitudinal collision ey energies. Due to our definition of the confining potential, the 
transverse excitation energies e±_ can take the possible values ej_ = e — ey = ^(2?7, + |m| + l) > 
of the discrete spectrum of the 2D oscillator |/xu;^p^. 

For the two-body interaction we choose a screened Coulomb potential 

V{r) = Vjj;e-'l'\ (6) 

already employed in ref. [211 I2S] for analyzing the ultracold scattering in cylindrical 
waveguides in the single-mode regime as = • The chosen potential ([6| depends on two 
parameters - the depth Vq < and the screening length ro > 0. Following the computational 
scheme already developed in ref. p6] , we implement different spectral structures of the atomic 
interaction by varying the single parameter Vq for a fixed length scale tq. Obvious advantages 
of the screened Coulomb potential ([6]) compared to the s-wave pseudopotential used in ref. 
[21] that is devoted to the multi-channel scattering of bosons are the following. First, by 
varying Vq one can vary the number of bound states of s-wave character in the interaction 
potential, and second, one can create new bound and resonant states of higher partial wave 
character. As a consequence we can consider not only bosonic but also fermionic as well as 
mixed collisions in a trap and including the case of higher energies. This will, as we shall see 
below, permit us to investigate new regimes and effects of multi-channel confined scattering. 
Performing the scale transformation 

r— >• — , e— >— , Vq ^ — and uj ^ — (7j 

ao Co Co uJq 

with the units oq = /i^//iVoro, eo = h? / iial, and ujq = e^/h, it is convenient to rewrite the 
equation ([s]) in the rescaled form 

-^V,2 + ^c.V + V(r)lv^(r)=6V^(r) , (8) 
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where V is now the correspondingly scaled potential and we fix /i = 1 and ro = 1 in the 
subsequent consideration. In the asymptotic region, \z\ oo, where the transverse trapping 
potential dominates the interaction potential, the axial and transverse motions decouple and 
the asymptotic wavefunction can be written as a product of the longitudinal and transverse 
0n,m(p, V^) components with 



/al)e*'"'^ (9) 

n\ ' 'a±- 
being the eigenfunctions of the transverse trapping Hamiltonian 

l,d^ 1 d 1 9\ 1 2 2 

with the corresponding eigenvalues e± = uj{2n + \m\ + 1) and the angular momen- 
tum projection m onto the z-axis. Here U^ix) are the generalized Laguerre polyno- 
mials and a_L = is the transverse oscillator length. The radial and azimuthal 
quantum numbers n and m of the 2D harmonic oscillator independently take the values 
n = 0, 1, 2, ...oo and m = 0, ±1, ±2, ... ± oo. Due to the axial symmetry of the transverse 
confinement and the spherical symmetry of the interatomic interaction ^ the angular mo- 
mentum component along the 2;-axis is conserved. Therefore, the problem ([s]) is reducible 
to a 2D one by separating the y9-variable and can be solved for every m independently. The 
quantum number n is a good one only in the asymptotic region — > oo and used for the 
definition of the initial (incident) asymptotic state i.e. channel 

C(r) = e^'=-0„,^(p,V.) , (11) 

of two infinitely separated atoms confined in a transverse state < pip\nm >. In addition 
to the quantum numbers n and m the asymptotic scattering state is defined also by the 
momentum /c„ of the channel 



2 \m\ 
kn = —\e-n - — , (12) 

which we express through the dimensionless energy e = e/{2uj) — 1/2. It is clear that the 
integer part of the dimensionless energy e coincides with the number He of open excited 
transverse channels. For n < rie we have {e — n — ^-^) > 0. The spherically symmetric 
interatomic interaction ^ mixes different transversal channels and leads to transitions n — >■ 
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n' between the open channels n, n' < ng, i.e. to transverse excitation/deexcitation processes 
during the coUisions. 

Assuming the system to be initially in the channel n, the asymptotic wavefunction takes 
at I -2 |— > +00 the form [29] 

rie 

^n,™(r) = e^'="^0„,™(p, ^) + X] [/nn' + ^^?^(^)/nn'] 6^'"' 0n',n. (p, (13) 

n'=0 

where /^^, and are the matrix elements of the inelastic scattering amplitudes for the 
even and odd partial waves, respectively, which describe transitions between the channels n 



and n' . For a bosonic (fermionic) collision just the symmetric (antisymmetric) part of (13) 
should be considered. 

It is clear that the scattering amplitude depends also on the index m which, however, 
remains unchanged during the collision due to the axial symmetry of the problem. Hereafter 
we consider only the case m = and the index is omitted in the following. 



Using the asymptotic wavefunction ( 13 ) and the total current conservation one obtains 



{Tnn' + Rnn' " K.n') = (14) 



n'=0 



for the inelastic transmission (reflection) coefficients T„„/ {Rnn')- kn {kn') is the initial (final) 



relative wave vector and n {n') the transverse excitation numbers according to eq.(12). We 
have [29] 

Tnn' = - n']^\6n,n' + fL' + fL'\' (15) 

"^n 

i?„„, = e[e-n']^|/l, + /„VP, (16) 

where G(x) is the Heavyside step-function. The transition probability Wnn', characterizing 
the transverse excitation/deexcitation, into a particular channel n' from the initial state n 
is given by the sum of the corresponding transmission and reflection coefficients 

Wnn' = Tnn' + Rnn'- (17) 

Due to the time reversal symmetry of the Hamiltonian we have Tnn' = Tn'n, Rnn' = Rn'n and 
Wnn' = Wn'n- The total transmission (reflection) coefficient T = J^n'^nn' {R = J2n' Rnn') 
is given by the sum of the transmission (reflection) coefficients of all the open channels. 



Eq.(14) leads to T + i? = 1. 
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III. NUMERICAL APPROACH 



To obtain the observable quantities T„„/ , and Wnn' of the scattering process, we 
have to calculate the matrix elements of the scattering amplitude / by matching the 



numerical solution of the Schrodinger equation ([8j) with the scattering asymptotics (13). 
To integrate this multi-channel scattering problem in two dimensions r and 6 [z and p) 
we adopt the discrete- variable method suggested in ref. |30j for solving nonseparable 2D 
scattering problems. This approach was applied in ref.[3T] to the case of a 3D anisotropic 
scattering problem of ultracold atoms in external laser fields. 

First, we discretize the 2D Schrodinger equation (8) on a 2D grid of angular {6'^}^]^ and 
radial {rj}^^^ variables. The angular grid points 9j are defined as the zeroes of the Legendre 
polynomial PNg{cos6) of the order Ng. Using the completeness property of the normalized 
Legendre polynomials which remains valid also on the chosen angular grid 

Ng-l 

J2 PlicOs9,)PiicOSej,)y^, = , (18) 

1=0 



where \j are the weights of the Gauss quadrature, we expand the solution of equation ([8 
in the basis fj{9) = J2f=o^ Pi{cos9)(P^^)ij according to 

Ng 



Hr,e) = lj^mu,ir) . (19) 
Here is the inverse of the Ng x Ng matrix P with the matrix elements defined as 



Pji = ^J\jPl{cos9j). Due to this definition one can use the completeness relation (18) in 
order to determine the matrix elements (P""'");^ as (P"''");^ = ^/\jPi{cos9j). It is clear from 



(19) that the unknown coefficients Uj{r) in the expansion are the values '?/'(r, 9j) of the two- 



dimensional wave function ipir, 9) at the grid points 9j multiplied by ^/\jr. Near the origin 



r ^ we have Uj{r) c^i r ^ due to the definition (19) and the demand for the probability 
distribution | ip{r,9j) p to be bounded. Substituting (19) into ([8| results a system of 
Ng Schrodinger-like coupled equations with respect to the A^e-dimensional unknown vector 
u(r) = {Af \,(r)}f« 

[H(°)(r) + 2(el- V(r))]u(r) = 0, (20) 

where 



V,y(r) = V{r,e,)5,,, = {V{r) + \uj^p]}6,,,, pj = rsm9, , (22) 



and I is the unit matrix. We solve the system of equations (20) on the quasi-uniform radial 
grid [32] 



r,=R-:^^ , J = 1,2,..., AT (23) 



^e^^i - 1 
- 1 

of grid points {rj} defined by mapping rj e (0, i? ^ +oo] onto the uniform grid Xj G (0, 1] 
with the equidistant distribution xj — = 1/A^. By varying and the parameter 7 > 
one can choose more adequate distributions of the grid points for specific interatomic and 
confining potentials. 



By mapping the initial variable r in Eq.(20) onto x we obtain 



[e(°)(x) + 2{el - V(r(x))}]u(r(x)) = , (24) 

with 



H\ 1 ^""^ 

^ ^ ^ 1=0 



where 



f{x) = . (26) 

The uniform grid with respect to x gives 6-order accuracy for applying a 7-point finite- 



difference approximation of the derivatives in the equation (24) . Thus, after the finite- 
difference approximation the initial 2D Schrodinger equation ^ is reduced to the system of 
N algebraic matrix equations 

3 3 

E ^^P^'J-P + + 2{el - V,}]u, + J2 ^+p^3+P = , J = 1, 2, AT - 3 

p=i p=i 

Uj + af^Uj^i + afuj_2 + + af^uj^^ = gj j = N -2,N -1,N (27) 

where each coefficient A^., is a A"e x Nq matrix, each a,- is a diagonal Nq x Nq matrix and 
each gj is a A'g-dimensional vector. Here the functions u_3, u_2, u_i and uq in the first 
three equations of the system (for j = 1,2 and 3) are eliminated by using the "left-side" 
boundary conditions: uq = and u_j = Uj {j = 1,2,3). The last three equations in this 
system for j = N, N — 1 and N — 2 are the "right-side" boundary conditions approximating 
at the edge points r7v-2, ''^at-i and r^v = -R of the radial grid, the scattering asymptotics 



(13) for the desired wave function u(rj). In order to construct the "right-side" boundary 
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conditions (27) at j = — 2, — 1 and we used an idea of ref.[3T] i.e. the asymptotic 



behaviour (13) at the edge points rN-2^^N-i and = R are considered as a system of 
vector equations with respect to the unknown vector /„„/ of the scattering amphtude for a 
fixed n. By ehminating the unknowns fnn' from this system we implement the "right-side" 



boundary conditions defined by Eqs.(27) at j = — 2, — 1 and A^ (see Appendix A). 
The reduction of the 2D muhi-channel scattering problem to the finite-difference bound- 



ary value problem (27) permits one to apply efficient computational methods. Here we use, 
in the spirit of the LtZ-decomposition [33], and the sweep method (or the Thomas al- 
gorithm [35J), a fast implicit matrix algorithm which is briefly described in Appendix B. 



The block-diagonal structure of the matrix of the coefficients in the system of equations(27) 
with the width of the diagonal band equal to 7 x A'e makes this computational scheme an 
efficient one. 



Solving the problem (27) for the defined initial vector fc„ and a fixed n from the possible 
set < n < rig we first calculate the vector function tlj{kn,r,9j). Then, by matching the 



calculated vector ip{kn, R, dj) with the asymptotic behaviour (13) at r = _R, we calculate the 
n-th row of the scattering amplitude matrix describing all possible transitions n ^ n' = 
0, 1, ...,ne. This procedure is repeated for the next n from < n < He- After calculating 
all the elements /„„/ of the scattering amplitude we obtain any desired scattering parameter 
T,R or W . 



IV. RESULTS AND DISCUSSION 



With the above-described method being implemented we have analyzed the two-body 
scattering under the transverse harmonic confinement for both cases of identical and distin- 
guishable colliding atoms. For confined scattering of identical atoms one has to distinguish 
the bosonic and fermionic cases. In the case of two colliding bosons the two-body wave 
function must be symmetric and only even scattering amplitude provides us with a nonzero 
contribution. First, we show that our result for the special case e < 1 of a single-channel 
scattering is in agreement at e ^ with the s-wave pseudopotential approach [13], and, 
particularly, reproduces s-wave CIR predicted and analyzed in refs. [131 [IH [211 EHl EH] 
for bosons. Then we extend our consideration to the multi-channel scattering £ > 1. We 
demonstrate that our results are in a good agreement in the limit of a long-wavelength trap 
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uo = 2ttc/X with the analytical expression given in [29] which has been obtained in 
the s-wave pseudopotential approach for the zero-energy limit. The range of validity of the 
analytical investigation in ref. PS] is explored. Next we present results for multi-channel 
scattering of two fermions under transverse harmonic confinement. For a fermionic colli- 
sion the two-body wave function is antisymmetric, i.e. only the odd scattering amplitude 
is nonzero. In the special case of single- channel scattering we reproduce the p-wave CIR 
for fermions [18j. These results are also in agreement with our previous investigations in 
refs. (2^ performed within a wave-packet propagation method |2S] • Finally we con- 

sider the confined multi-channel scattering of two distinguishable atoms. In this case both 
even and odd amplitudes contribute to the scattering process. 




FIG. 1: s- and p-wave scattering parameters (red) and Vp (black) as a function of the potential 
depth Vo for the free-space effective potential V{r) + l{l + l)/(2r^). Divergences correspond to the 
appearance of new bound (s-wave scattering) or shape resonant (p-wave scattering) states in the 
effective potentiaL All quantities are given in the units ([7|. 

For modeling different interatomic interactions in the subsequent sections we vary the 
depth Vo of the potential (|6]) in the wide range — 10 < Vq < for a fixed width tq = 1. In 
free space, this potential being superimposed with the centrifugal term /(/ -|- l)/(2r^) makes 
an effective potential which may support some even or odd bound states depending on the 
value of the quantum number I and the parameter Vq. For / ^ there might be also some 
shape resonances for certain relative energies. These shape resonances may enhance strongly 
the contribution of Z 7^ partial waves in the energy domain where one would have expected 
a pure / = scattering. 
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It is known that in the zero-energy hmit, when the scattering process does not depend on 
the details of the potential, the collision can be described by a single parameter: the s-wave 
scattering length = —limk-,Q tan 6s{k)/k for / = (bosonic collision) and the p-wave 
scattering volume Vp = —limk^o tan 6p{k)/k^ for / = 1 (fermionic collision). For sufficiently 
low collision energy the contribution of the partial waves with larger I can be neglected. 
In Figjl] we have plotted a<j and Vp in the region —10 < Vq < 0. Figjl] demonstrates the 
rich spectral structure of the chosen form of the interatomic interaction: the scattering 
parameters and Vp can be positive or negative and they diverge for the values of Vq 
corresponding to the appearance of new bound states. In the case of p-wave scattering the 
increase of the depth Vq of the potential first leads to a shape resonance which approaches 
zero energy and finally transforms to a p-wave bound state. 



A. Mult i- channel scattering of bosons 

For two bosons colliding in a transverse harmonic confinement the scattering wave func- 



tion is symmetric with respect to the exchange z —z, i.e. /°^, = in Eq.(13). We 
consider multi-channel scattering with the dimensionless energy e = e/{2uj) — 1/2 < 4 (12) 
permitting the collisional transverse excitations/deexcitations n — > n' < 4 up to four open 
channels (four-mode regime). For comparison with analytical results [HI [21] obtained in the 
s-wave pseudopotential approach, we have extracted the effective quasi-lD coupling constant 



giD = limk^o Re{fQQ{k)} / Im{fQQ{k)}k / jj, as well as the transmission coefficient T = Too(15) 
and the scattering amplitude /qq([I3|), as a function of the scattering length as in the single- 
mode regime (0 < £ < 1). The calculated parameters are presented along with the analytical 
results in Fig|2] for u = 0.002 and the longitudinal relative energy e\\ = e — e± = 0.0002. 
Fig|2]^a) shows the coupling constant gir, as a function of the scattering length a^. Our 
numerical result clearly exhibit a singularity at as/a± ~ 1/C with C = —^(1/2) = 1.4603.., 
which corresponds to the well-known s-wave CIR [13[[T1]. In Fig|2|^b) we present the trans- 
mission coefficient T versus a^. The transmission coefficient T goes to unity (total trans- 
mission) when as tends to zero (i.e. no interaction between the atoms), while at the CIR 
position, it exhibits the well-known minimum (blocking of the atomic current by the CIR). 
Fig|2|c) shows the scattering amplitude /qq as a function of a^. The amplitude approaches 
zero at as/a_L = and —1 at the CIR position , which results in total transmission and total 
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reflection respectively. In general the presented values of /qq are in very good agreement 
with the analytical results. 




-3 -2-10123 



FIG. 2: (a) The effective quasi- ID coupling constant gn), (b) the transmission coefficient T and (c) 
the scattering amplitude /qq as a function of the scattering length Ug for one-channel scattering of 
two bosons under a harmonic confining potential with uj = 0.002 for a longitudinal relative energy 
ey = 0.0002 together with the analytical results obtained for the s-wave pseudo-potential zero- 
energy limit(solid curves ) [131 [HI [29] . The constant Vq is varied in the region —2.30 < Vq < —0.32. 

In the multi-mode regime our results are in good agreement with the analytical ones 
obtained within the s-wave pseudopotential approach in the zero-energy limit [22] for long- 
wavelength traps (tight confinement: to = 27rc/A —^0). In Figjsjwe present our scattering 
amplitude /qq along with the analytical results as a function of the dimensionless energy 
e foT UJ = 0.0002 and Vq = —3.0. Apart from energies close to the channel thresholds, 
the real and imaginary part of the scattering amplitude /qq show a monotonous behaviour: 
i?e(/Qo) is monotonically increasing and approaching zero asymptotically whereas /m(/oo) 
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decays monotonically and also approaches zero for large values of e. The peak structure 
located at integer values of e is due to the resonant scattering once a new previously closed 
channel opens with increasing energy. There is a good agreement between our results and the 
analytical ones given by Eq.(6.9) in ref. [29] for the complete range < e < 4. However, we 
encounter major deviations with increasing uj except for narrow regions close to the channel 
thresholds (as e approaches integer values) for the real parts of the scattering amplitudes 
(see Figji]). These deviations are most presumably due to the energy dependence of the 
s-wave scattering length which is neglected in ref. |29] . 




FIG. 3: The scattering amphtude /qq (dots) as a function of the dimensionless energy e along with 
the analytical results (solid curves) for uj = 0.0002 and Vq = —3.0 (a^ = —8.95). 

For the multi-channel regime {e > 1) we find a strong dependence of the total transmission 
coefficient T on the population of the initial state n, except for the case cts/aj, ^ of 
noninteracting bosons in free space. Figjs] shows the calculated transmission coefficients(15) 
as a function of as/a± for uj = 0.002 for the (a) two-mode regime with e = 1.05 and (b) the 
three-mode regime with e = 2.05. Note, that all these cases correspond to near-threshold 
collision energies if the maximal integer is subracted from e. Similar to the single-mode 
regime, T goes to unity (total transmission) when tends to zero. We also encounter a 
minimum. However, the value of the transmission at the minimum is not zero anymore in 
the multi-mode regime, the larger the number of open channels, the position of the minima 
will be more shifted to the left. For a fixed value of the ratio as/a± 7^ 0, a lower initially 
populated transverse level n, leads to a larger total transmission. 

In Fig|6] we present the transmission as a function of the dimensionless energy e for the 
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FIG. 4: The scattering amplitude /qq (dots) as a function of the dimensionless energy e along with 
the analytical results (solid curves) for lo = 0.02 (red), uj = 0.04 (blue) and u = 0.06 (black). The 
amplitudes have been calculated for Vq = —3.0 (a^ = —8.95). 
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FIG. 5: The total transmission coefficient T for a bosonic collision in a harmonic confinement 
u; = 0.002 as a function of as/a± for (a) two-open channels with e = 1.05 being initially in the 
transverse ground state n = (red) and first excited state n = 1 (black), and (b) three-open 
channels with e = 2.05 being initially in the ground state n = (red), first excited state n = 1 
(black) and second excited state n = 2 (blue). 

two cases of the system being initially in the ground n = (a) or first excited n = 1 {h) 
transverse states, for different values of the ratio as/a±. With increasing collision energy 
the transmission coefficient exhibits a nonmonotonic behaviour: we observe a sequence of 
minima and peaks. For < e < 1 the value of T at the minimum is zero. As the number 
of open channels increases with increasing collision energy e, the transmission values at the 
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corresponding minima increase strongly. The corresponding transmission peaks T = 1 are 
located at the channel thresholds. The shape of the transmission Valleys' in between two 
integer values of e as well as the positions of the minima strongly depend on the ratio cis/a_L 
which can be changed by varying the strength of the interatomic interaction Vq or the trap 
frequency uj. Both, changing Vq and/or uj leads to energetical shifts of the bound states (in 
particular for the excited transversal channels) of the atoms in the presence of the confine- 
ment. If the collision energy coincides with a bound state of the corresponding closed channel 
we encounter an occupation of the closed channel in the course of the scattering process, 
i.e. a Feshbach-resonance occurs. The interpretation of the minimum of the transmission T 
in terms of a Feshbach-resonance at the point as/ai. = 1/C for the zero-energy limit was 
provided in refs. [I3l [H] , where it was shown that the origin of the CIR is an intermediate 
occupation of a bound state belonging to an excited transverse (closed) channel. 

To demonstrate that the above-discussed behaviour (minima) of the transmission co- 
efficient T{e) in certain regions of e is due to Feshbach resonances we have analyzed the 
probability density of the scattering wave function of the atoms in the trap. In the single- 
mode regime and in the zero-energy limit we encounter the well-known CIR: Fig j7] shows the 
corresponding probability density |'?/'(x, z)p for an initial transverse ground state n = and 
e = 0.05 as well as as/a± = 0.68. For small \z\ one observes additional two pronounced peaks 
along the transverse (x-) direction corresponding to the occupation of the bound state (with 
the binding energy e^^^ ~ 0) in the first excited closed channel of the transverse potential. 
The probability density tends to zero as 2; ^ +00. This leads to a zero of the transmission 
T{e) for £ ^ (see Figj6|^a)) corresponding to the zero-energy CIR. 

In Figjsjwe show the probability densities \iIj{x,z)\'^ at as/a± = -1-4.39 for several values 
of the dimensionless energy e for (a) the single- and (b) two-mode regimes for collisions 
with initial transverse state n = 0. The corresponding probability density exhibits for small 
values of 1 2; I additional two (for the single-mode regime) and four (for the two- mode regime) 
pronounced peaks with respect to the transverse (x-) direction as e approaches the CIR- 
position. This demonstrates the occupation of bound states of higher, namely first excited 
(with binding energy e^^i) and second excited (with binding energy £^=2) channels in the 
course of the scattering process. The corresponding transmission values are also indicated in 
Figjs] Zero transmission is observed also for as/a± = -1-4.39 and e = 0.75 with no probability 
density being present for large positive values of z, see Figjsj^a). 
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FIG. 6: The total transmission coefficients T for bosonic collisions as a function of the dimensionless 
energy e for the two cases of the system being initially in the ground n = (a) and first excited 
n = 1 (b) transverse states, for several ratios of as/a± and uj = 0.002. The black curve corresponds 
to as/a± = 1/C for which the zero-energy CIR in the single-mode regime is encountered. 



T=0.0 




FIG. 7: The probability density |V'(x, z)\'^ for bosonic collisions as a function of x and z at as/a± = 
0.68 for e = 0.05 - zero-energy CIR . The corresponding transmission values are also indicated. 
The result has been obtained for uj = 0.002. 

The energy of the bound state that leads to a minimal transmission due to a 

resonant scattering process, changes with varying as/a± as follows. It is below e = for 
< as/a± < 0.68 and consequently no minimum is encountered for T{e) in the range 
< e < 1 (see fig. [6]). At the position of the zero-energy CIR it is located just above the 
threshold e = leading to a zero transmission for e ^ 0. For as/a± > 0.68 the bound 
state energy is somewhere in between the channel thresholds e = and e = 1 whereas for 
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FIG. 8: The probability density \tp{x, z)p for bosonic collisions as a function of x and z at as/a± = 
+4.39 for two cases of the single-mode regime (a) and two-mode regime (b) with different values 
of e. The corresponding transmission values are also indicated. All subfigures are for u = 0.002 
and 77, = 0. 
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CLs/cL± < it is below but close to e = 1 leading again to a corresponding minimum of 
T. For both cases, as/a± < and as/a± > 0.68, an increase of as/a± leads to a narrow 
transmission well and the corresponding transmission minimum is shifted towards the next 
higher channel threshold. The dependence of e^^i on the parameter as/a± is in agreement 
with the pseudopotential analysis given in ref. |14j (see Fig. 2 of this ref.). 

In Figj9]^a) we show the transition probabilities Wnn' as a function of as/a± for to = 0.002 
and 6 = 3.05 corresponding to four open channels. We observe that the probability of 
remaining at the same initial state, Wnn (i-e. elastic scattering) is in the complete range 
of the ratio as/a± much larger than the probability of a transition into a different state, 
Wnn' (i-e. inelastic scattering). With increasing n or n' the inelastic transition probabilities 
Wnn' increase but the elastic probabilities Wnn decrease. In an inelastic (elastic) collision 
Wnn' (Wnn) goes to zcro (uuity) as as tends to zero. Wnn' (Wnn) possess a maximum 
(minimum) at the resonance position as/aj_ ^ 0.35 consistent with the minimum of T{e) 
at e = 3.05. It is instructive to see how the distribution of the initial flux among the open 
channels changes due to pair collisions as a function of the collision energy. Fig|9]^b) shows 
the transition probabilities Wnn' as a function of the dimensionless energy up to four open 
channels. The probability of elastic scattering remains larger than that of inelastic scattering 
in the complete range of the energy. For two open channels the elastic collision probability 
Wnn is independent of the initial state (Woo = Wu). For a higher number of open channels 
Wnn is decreasing with increasing initial value of n. Near the thresholds, the probabilities 
of the inelastic (elastic) transitions Wnn' (Wnn) go to zero (unity). 



B. Mult i- channel scattering of fermions 

In this section we focus on fermionic collisions in harmonic traps. In this case the 
interatomic wave function is anti-symmetric with respect to the interchange of the two 



fermions and the even amplitude f^^, in eqs.(13) is zero. We have analyzed the multi- 
channel scattering of fermions up to four open transverse channels for different interatomic 
interactions, by varying the potential strength Vq in the vicinity of the value Vq = —4.54 



(see Fig.l) generating a resonant p-wave state in free space. Fig. 10 shows correspond- 
ing results for the single-mode regime. Fig. [lO|^a) shows the mapped coupling constant 
911? — li™A:-»o Ii^{f°{k)}/Re{f°{k)}jj,k [HI [26] as a function of the p-wave scattering length 
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FIG. 9: (a) The transition probabilities Wnn' as a function of as/a± in the four mode regime 
for e = 3.05. (b) The calculated transition probabilities Wnn' as a function of e up to four open 
channels for as/a± = —30.08. uj = 0.002 for both subfigures. 



ttp = yVp- The mapped coupling constant g^^y goes to zero at the position of the mapped 
CIR llSj, which for uj = 0.002 and ey = 0.0002 is equal to ap/a± = —0.31. The position of 



the mapped CIR obviously depends on the values of uj and e. In Fig 10 b) we have plotted 



the transmission coefficient T as a function of ap. The transmission coefficient exhibits a 
minimum and an accompanying well (i.e. the blocking due to the resonance) at the position 
of the CIR, and tends to unity (i.e. total transmission) for ap far from the CIR-position. For 
larger u the well becomes wider, and its minimum is shifted to the left, see also ref. [24]. Fig. 
[Toj^c) shows the scattering amplitude /qq as a function of ap for to = 0.002 and e\\ = 0.0002. 
The amplitude approaches zero far from the CIR position and —1 at the CIR position , 
which results in total transmission, T = |1 + /qq^ — 1 (i.e. the fermions do not scatter each 
other) and total reflection, T = |1 + /qoP (i.e. strongly interacting and impenetrable 
fermions), respectively. 



In Fig, 11 we present the total transmission coefficient as a function of ap in the two- mode 



regime for e = 1.05 (a) and three-mode regime for e = 2.05 (b) for u = 0.002. Similar to the 
single-mode regime, T exhibits a minimum and accompanying well, however the position of 
the minimum is shifted to the right and the value at the minimum is nonzero. For the lower 
degree of transversal excitation, we observe a deeper and narrower transmission well. With 
increasing energy the transmission well becomes wider and more shallow and its position is 



shifted to larger values of ap/ a±. This is also demonstrated in Fig, 12 , where the transmission 
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FIG. 10: (a-b) The mapped coupling constant g^^^ and the transmission coefficient T as a function 
of the scattering length ap for single-channel scattering of two fermions under a harmonic confine- 
ment Lo = 0.002 and for the longitudinal energy = 0.0002 (red), u = 0.02, = 0.002 (black) 
and uj = 0.06, ey = 0.002 (blue), (c) The scattering amplitude /qq as a function of the scattering 
length Up for uj = 0.002 and e|| = 0.0002. The corresponding constant Vq is varied in the region 
-4.54 <Vo< -4.47. 
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coefficient is plotted as a function of e for several values of a^/ a±_ for the two cases of being 
initially in the ground n = (a) and the first excited = 1 (b) transverse states for 
uo = 0.002. We observe that for any number of open channels, T exhibits a minimum for 
some value of ap/a±. With increasing ap/a± the transmission well becomes more shallow 
(except for the single-mode regime) and wider. In contrast to the bosonic case, there is no 
specific threshold behaviour. This is a consequence of the fact that the relative motion does 
not feel the interatomic interaction in the closed channels which are strongly screened by 
the centrifugal repulsion playing a dominant role for near-threshold collision energies or in 
other words: We encounter a weak coupling of the different scattering channels. 
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FIG. 11: Total transmission coefficient T in a fermionic collision as a function of Up for uj = 0.002 
for (a) two-open channels with e = 1.05 being initially in the ground state n = (black) and first 
excited state n = 1 (red), and (b) three-open channel with e = 2.05 being initially in the ground 
state n = (blue), first excited state n = 1 (black) and second excited state n = 2 (red). 



In Fig, 13 we show the transition probabilities Wnn' as a function of a^/ a± for four open 
channels. We see that the probability of an elastic scattering process (i.e. to remain in the 
same transversal state) is much larger than that of an inelastic collision (i.e. the transition 
to a different transversal state). For an elastic collision the probability Wnn shows a mini- 
mum and corresponding well which becomes wider and more shallow with increasing initial 
quantum number n (i.e., with the population of a higher excited initial transversal state). 
For an inelastic collision Wnn' (^ 7^ n') exhibits a peak which becomes less pronounced as 
the quantum numbers n or n' increase. 
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FIG. 12: Total transmission coefficient T for a fermionic collision as a function of the dimensionless 
energy e for the two cases of the system being initially in the ground n = (a) and first excited 
n = 1 (b) transverse states, for several ratios of a,p/a±. We have used uj = 0.002. 



FIG. 13: Transition probabilities Wnn' in fermionic collisions as a function of ap/a±_ for four open 
channels, uj = 0.002 and e = 0.0142 are employed. 

C. Mult i- channel scattering of distinguishable particles 

In this section we analyze the multi-channel scattering of two distinguishable particles 
in a harmonic trap with the same trap frequency uji = uj2 = oj allowing the separation of 
the cm. motion. Such a case corresponds to the different atomic species confined by the 
same potential as it is the case e.g. for two different isotopes in the same optical dipole 
trap. The two-body scattering wave-function of the distinguishable atoms does not possess 
a well-defined symmetry with respect to the reflection z — > —z, i.e. both s- and p-state 
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contributions (/*^ and f°) must be taken into account for the scattering amplitude. Fig, 14 
shows our results for the transmission coefficient T in the single-mode regime (0 < £ < 1), 
which is plotted as a function of the tuning parameter — Vq of the interparticle interaction. 
In general the scattering process can not be described by a single scattering length as or 
for this case. In regions with a negligible p-wave contribution (see the regions in Figjljwith 
Vp ^ 0), T exhibits a behavior similar to bosonic scattering. We observe the well-known 
s-wave CIRs which lead to zeros of the transmission T at the positions as/a± = 1/C and 
tend to unity when goes to zero (together with Op). In regions, where and Op are 
comparable we observe the effect reported in ref. pi]: remarkable peaks of the transmission 
T = |1 + /oo + /ooP ~^ l i-^- almost complete transmission in spite of the strong interatomic 
interaction in free space. This is the so-called dual CIR: Quantum suppression of scattering 
in the presence of confinement due to destructive interference of odd and even scattering 
amplitudes. Equally minima of T = |1 + /oo + /ooP ~^ interference of even 

and odd scattering amplitudes under the action of the transverse confinement can occur 
(see Fig, 14). Complete transmission corresponds to foo + foo = —^ — 'i while total reflection 
corresponds to /qq + /oo — Fig' 15 shows the corresponding amplitudes /qq + /qo as a 
function of —Vq. 




FIG. 14: Transmission T as a function of the depth —Vq of the potential Q for two distinguishable 
atomic species for uj = 0.002 for en = 0.0002. Units according to eq.(l7l). 



Fig 16 shows the total transmission coefficient versus — Vq in the three-mode regime for 
UJ = 0.002 and e = 2.05. Similar to the single-mode regime, when ap is negligible compared 
to a^, T behaves analogously to the case of a bosonic collision, and tends to unity (complete 
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FIG. 15: Scattering amplitude as a function of the depth —Vq of the potential ([6| for two distin- 
guishable atomic species for u> = 0.002 and e|| = 0.0002. Units according to eq.([7|. 



transmission) when tends to zero, while at the s-wave CIR position, it exhibits a minimum 
with a nonzero value. For a lower degree of transversal excitation of the initial state, we 
encounter a larger transmission coefficient. For the same reasons as in the single- mode 
regime, we observe in the regions of Vq where the p-wave scattering length is comparable 
to ttg, sharp peaks of T. However in contrast to the single-mode regime we do not observe 



complete transmission i.e. T ^ 1 (see e.g. T at Vq = —8.85 in Fig, 16) 




FIG. 16: Transmission coefficient T as a function of the depth — Vq of the potential ^ for two 
distinguishable atomic species for oj = 0.002 and e = 2.05 in the three-mode regime, being initially 
in transverse ground state (blue), first excited state (red) and second excited state (black). Units 
according to eq.Q. 
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Fig 17 presents the transmission coefficients as a function of the dimensionless energy e 



for uj = 0.002. Fig, 17 a) shows the results for several ratios of as/a± if \ap\ is small compare 
to \as\. Here the system is initially in the transversal ground state. The behaviour of the 
transmission coefficient is similiar to the case of two-boson scattering (see|6|-a). FigjlT^b) 
shows the transmission for the scattering of two distinguishable atoms together with the 
results for bosonic- (blue) and fermionic- (black) collisions, for — Vq = —4.505 when both 
scattering length and are large. In the limit e the behaviour of the transmission for 
the distinguishable and the bosonic case are very similar while in the vicinity of the energy 
of the p-wave shape- resonance e ~ 0.7, apart from a small shift to larger energies, we can 
find a complete coincidence of the transmission behaviour for the case of distinguishable and 
fermionic scattering. At the threshold energy e = 1 the latter two transmission curves cross. 
Figjl7|(c) shows the transmission coefficient versus e for — Vq = —4.505 initially occupying 
the ground state n = (blue), the first excited state n = 1 (black), the second excited state 
n = 2 (red) and the third excited state n = 3 (green). In the limit — > 0, T drops rapidly 
to zero. Apart from the single mode regime, it is for >> approximately constant with 



In figjlSjwe present the transition probabilities Wnn' as a function of — Vq for four open 
channels. For the regions of Vq where is negligible compared to a^, a very good agreement 
of the transmission behaviour for distinguishable and bosonic atoms is observed. In regions 
where and Op are comparable, similar to the fermionic case, we observe narrow and deep 
wells (for elastic collision) and narrow as well as strongly pronounced peaks (for inelastic 
collision) . 



V. SUMMARY AND CONCLUSIONS 



We have analyzed atomic multi-channel scattering in a 2D harmonic confinement. Iden- 
tical bosonic and fermionic scattering as well as scattering of distinguishable atoms in traps 
with the same frequency for the different species have been explored. Equal frequencies 
allows us to separate the cm. and relative motion. 

Firstly we reproduced the well-known s-wave CIR for bosonic collision in the single-mode 
regime [131 El Ell ESI EEl ES]. Next bosonic collisions in the multi-channel regime including 
elastic and inelastic processes i.e. transverse excitations and deexcitations have been in- 
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FIG. 17: Transmission coefficient T as a function of the dimensionless energy e for scattering of 
two distinguishable particles in a harmonic confinement a; = 0.002 (a) for several ratios as/a± in 
the case where Op is negligible compared to Ug, (b) for — Vb = —4.505 when both scattering lengths 
Os and Op are comparable (here the transmission coefficients for bosonic- (blue) and fermionic- 
(black) collision are also provided) and (c) for —Vq = —4.505 when the system is initially in the 
ground transversal state (blue), in the ffi'st excited state (black), in the second excited state (red) 
and in the third excited state (green). 
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FIG. 18: Transition probabilities Wnn' for quasi-lD scattering of distinguishable particles in a 
harmonic trap uj = 0.002, as a function of — Vq for four open channels and e = 2.05. Vq is given in 
units of eq.Q. 

vestigated. Transmission coefficients as well as transition probabilities for energies covering 
up to four open channels are reported on. It is shown that the transmission coefficient as 
a function of the scattering length exhibits a minimum at the CIR position. Except for 
< as/a± < 1/C in the ffist channel, the transmission curves show, with varying energy, a 
minimum and accompanying well. For a single open channel the value of the minimum of 
the transmission is zero. Increasing the degree of transverse excitation leads to an increase 
of the transmission minima. A distinct threshold behaviour is observed at every e = n. 

For as/aj_ < the accompanying transmission well is more pronounced compared to the 
case as/a± > and it is also deeper in the case of several open channels. The position of 
the minimum is closer to the upper threshold again when compared to the case as/a± > 
0. With increasing ratio as/a± for either as/a± > or as/a±_ < 0, the corresponding 
transmission well becomes also more pronounced i.e. we encounter a deeper and narrower 
well for more open channels and the energetical position of the minimum moves to the 
neighboring upper threshold. Our results for the transition probabilities Wnn' show that the 
probability of remaining in the same state (elastic collision) is larger than transitions into 
a different transversal state (inelastic collision). In the case of elastic collisions (inelastic 
collisions) it has to go to unity (zero) as as tends to zero. The inelastic transition probability 
Wnn' increases as the numbers n or n' increase while the elastic transition probability Wnn' 
decreases. 
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Our next focus has been the multi-channel scattering of fermions addressing up to four 
open channels in the waveguide. In the regime of a single open channel we reproduced the 
well-known p-wave mapped CIR in the zero-energy limit [TSl (SHj. We have analyzed the p- 
wave CIR dependence on the collision energy e and the trapping frequency uj. Consequently 
the multi-mode regime has been explored by determining the behaviour of the transmission 
as well as the transition probabilities. The dependence of the transmission on the p-wave 
scattering length exhibits a minimum with zero transmission and an accompanying well for 
the single-mode regime. Increasing the energy, this well becomes more shallow (except for 
the single-mode regime), wider and its position is shifted to the region of higher energies. 
In the multi-mode regime the transmission well depends also on the degree of the initial 
transverse excitation. Increasing the degree of initial transverse excitation the transmission 
well becomes shallower and wider and the transmission is overall increased. For a fixed 
number of open channels the transmission exhibits as a function of the energy a minimum 
for some value of the p-wave scattering length. Increasing a^ja^ leads to a shift of the 
transmission minimum to higher energies and an increasingly shallower well. This holds 
except for the single mode regime in which the minimal value is zero. 

In contrast to the bosonic case, we do not observe a distinct threshold behavior. The 
transition probabilities Wnn' for fermionic collisions show that elastic collisions are more 
probable than inelastic ones. The probability Wnn for an elastic process shows a well which 
becomes wider and shallower as the initial population of the excited states n increases, while 
for an inelastic collision it exhibits a peak which becomes smaller as the channel numbers n or 
v! increase. By varying the potential parameter Vq in the viccinity of the value Vq = —4.54, 
we found a set of shape-resonances. 

Finally we analyzed the multi-channel-lD scattering of distinguishable atoms in the trap. 
Here both s- and p-wave contributions have been taken into account. We have studied the 
transmission coefficient in the single-mode regime and reproduced the CIR [131 EH [ISl EH] 
as well as the dual-CIR [2ll ESI EH] corresponding to total reflection and transmission, 
respectively. In the multi-mode regime the transmission versus — Vq shows except for the 
regions where is comparable with a behaviour analogous to the situation of bosonic 
scattering. We encounter T ^ for 0. At the position of the CIR, T exhibits a 

minimum with a nonzero value. The lower is the transverse excitation of the initial state, 
the larger is the transmission coefficient. In regions of Vq where the p-wave scattering length 
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Op is comparable to Os, T exhibits sharp peaks or dips near the dual-CIRs due to comparable 
contributions of both s- and p-waves. In contrast to the single-mode regime there is no a 
complete transmission at the dual CIR points. Our results for the transmission coefficient 
versus energy show that in the limit of zero longitudinal relative energy the distinguishable 
atoms behave like bosons, while in the vicinity of the position of a shape-resonance, they 
behave like fermions. For larger energies e >> 1 T is close to unity. Finally we have analyzed 
the transition probability as a function of — Vq. For values of Vq where is negligible 
compared to we find an excellent agreement with the results for bosonic collisions. In 
regions where these two scattering parameters are comparable, we observe, similar to the 
fermionic case, sharp downward- (for elastic collision) and upward- (for inelastic collision) 
peaks but with different values. We conclude with the general statement that our multi- 
channel scattering results in waveguides are of immediate relevance to cold or ultracold 
atomic collisions in atomic waveguides or impurity scattering in quantum wires. 
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VII. APPENDIX A. BOUNDARY CONDITIONS 

In the limit r tends to zero, u goes to zero. This boundary condition can be satisfied 
easily by putting u(r = 0) = and u j = uj in the vicinity of r = 0. The latter is needed 
for approximating the derivatives ^■u(r(x)) and ^■u(r(x)) near the point r = 0. 



The boundary condition for u approximating the asymptotic form (13) at large r can be 
written in the form 

u, + af^Uj^i + afuj_2 + af ^u,_3 + afu^^^ = g, j = iV - 2, iV - 1 , AT (28) 

where a^s are diagonal Nq x A^^^ matrices and gj is a A^^e-dimensional vector. The above 
boundary-conditions are constructed by eliminating the unknown amplitudes from the 
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asymptotic equations (13) written for a few rjs neighbouring to the point = R. This 



gives the values for the coefficients a^s and gj. In general for up to four open channels we 
have 



a 



(0 



and 



= (-1)^ (T^'^^rn^^wrA 1^ ^ML^^A\-^T\~\-^.,\) ^ / = 1, 2, 3, 4 (29) 



(30) 



0n (pT /) 

For the case of bosonic (fermionic) collisions we must consider just the even (odd) part of 
the equation (30), i.e. we just need to replace the term e^^"^]-i by cos knzj^i {i sin knz"!_i). 



Here 0n(p) = 4>n,o{p,v) is the eigenfunction of the transverse trapping Hamiltonian (see 
Eq.([9|), = Tj sinOm, = rj cosOm, n is the channel number of the initial state, and Ug 
is the number of transversely excited open channels. The complex numbers Tj'^^s are given 
through 

Tj'™ = Si^o + Si,i{af + + cf) + (5z,2(a767_ 1 + afc]l, + bfcf.,) + 5i,,a^bf_,c';'_^ (31) 
If ne < 1 (a single-mode regime) a™ = 0, otherwise 

- l-0,i - U + + Cj ''-Oj-l + \0j + Cj + Cj_ij^o,i-2 - Cj_i"0,j-3/ X 

if ne < 2, 6™- = 0, otherwise 

6f = {H™. - (1 + cf )H- _i + cf .J X {S™. _i - (1 + cf )S- _2 + H™. _3}-' (33) 



if ne < 3, c™ = 0, otherwise 



-1 



(34) 



there SJ^^- 



' 'JZ(P^' ^^'^ ^""^ ^^"^^^ eq.(|12|) 



VIII. APPENDIX B. FAST IMPLICIT MATRIX ALGORITHM 



Following the idea of the LU-decomposition [53] and the sweep method [31] (or the 
Thomas algorithm ^35]), we search the solution of the system of vector equations with 



any coefficient being a. Nq x Ng matrix (27) in the form 

„ ^r^^\. 



j = l,...,iV-3 



(35) 
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Here we define u{r{xj)) as Uj for simplicity. Tlie CjS are unknown Ng x Ng matrices. To find 
the solution, first we should calculate the unknown Cj matrices. The plan is the following: 



Due to eq.(35) we have 



u 



r'(^) 11 -I- r^-"' 11 -I- r^-^> n 



^(2) 



(3) 



J-p 



j-p+35 



p = l,2,3 



(36) 



then one obtains 



and 



„ _ M^) „ I r^(2) , ^(3) 



U,-2 



(2) 

i-2 



/^(l) /^(2) I ^voj ^vi; ^yo) 



.(3) 



(1) r-(3) 



(37) 

(38) 



u 



j-3 



r^(i) ^(1) 



^(1) ^(2) 



/^(2) ^(1) I /^KOj 
S-2S-I + S-2 



(3) 



+ 



^(1) ^(1) ^(2) 



+ 



^(1) ^(3) 
■^(1) ^(1) ^(3) 



^(2) ^(2) 
^(2) ^(3) ■ 



U 



J+1 



u 



J+2 



(39) 



By substituting Uj defined by Eqs.(37 39) into Eq.(27) one can calculate Uj in terms of 
Uj+i, Uj4.2 and Uj+3. Then, by comparing with Eq.(35) we find a recurrence formula for 
calculating the unknown matrices Cj: 



DC) 



(1) 



^i-3 



+ 



■^(1) /^(l) /^(2) , ^(1) ^(3) ^(2) ^(2) ■ 

^,-^cfl^fl.+cfl^+^_^fl,+^,, 



(40) 



-DC 



(2) 



^i-3 



^(1) ^(1) ^(3) 



r^(2) ^(3) 



^(1) ^(3) 
-2<-'i-2<-'i-l 



^j+2 



(41) 



and 



DC 



(3) 



(42) 



Here 



^i-3 



^(1) ^(1) ^(1) 



^(1) ^(2) ^(2) ^(1) ^(3) 
S-3S-2 + S-3S-I + S-3 



■J-2 



Cf^^f}, + C^, +Aj_,C;;^,+A^ + 2(e/-l^) 



.(2) 



(1) 



(43) 



By using the left-side boundary conditions and Eq.(35) one can calculate the Cj matrices for 
j = 1, 2 and 3. Then by using Eqs.(43 40) we calculate all the matrices Cj. Subsequently by 
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using the right-side boundary conditions ( 28 ) and recurrence formula ( 35 ) we first calculate 
Uj foT j = N — 2,N — 1 and and then for j = 1, ...,N — 3. 
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